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We derive an expression that allows for the unambiguous evaluation of the overlap between two 
arbitrary quasiparticle vacua, including its sign. Our expression is based on the Pfaffian of a skew- 
symmetric matrix, extending the formula recently proposed by [L. M. Robledo, Phys. Rev. C 79, 
021302(R) (2009)] to the most general case, including the one of the overlap between two different 
blocked n-quasiparticle states for either even or odd systems. The powerfulness of the method is 
illustrated for a few typical matrix elements that appear in realistic angular-momentum-restored 
Generator-Coordinate Method calculations when breaking time-reversal invariance and using the 
full model space of occupied single-particle states. 



I. INTRODUCTION 

The evaluation of kernels in projection and more gen- 
eral applications of the Generator Coordinate Method 
(GCM) based on quasiparticle vacua requires the calcu- 
lation of the overlap between two different quasiparticle 
vacua. Its evaluation presents a long-standing technical 
challenge: the standard expression for this overlap, the 
so-called " Onishi formula" p] provides the square of the 
(complex) overlap only. As a consequence, the overall 
sign of the overlap is not determined, or, equivalently, its 
phase is determined up to integer multiples of ir only. 

One possible solution to the problem was proposed by 
Neergard and Wiist [3J, but the practical application of 
their technique becomes cumbersome in realistic appli- 
cations and has been rarely used in practice. Notable 
exceptions are Refs. [2111]. Many groups have resided to 
determine the phase through a kind of Taylor expansion 
that allows "to follow the overlap" when the kernel can 
be connected in small steps to a known reference over- 
lap (see, for example, Refs. [51 [5]). In practice this can 
become very cumbersome or even impossible when the 
phase is rapidly changing or when there is no symmetry 
that establishes a reference phase. 

It was pointed out by Robledo in Ref. [7] that tech- 
niques for the manipulation of matrix elements between 
Fermionic coherent states that are well-known in field 
theory allow to express the overlap between two quasipar- 
ticle vacua, including its sign, as the so-called Pfaffian of 
a skew-symmetric matrix. 1 In a more recent paper, Rob- 
ledo [5] has also worked out the practical implementation 
of this idea for the unambiguous evaluation of the over- 
lap between fully-paired quasiparticle vacua, including 
the limit where some of the "pairs" consist of fully occu- 
pied single-particle states. Even more recently, Bertsch 
and Robledo [TU] also investigated the case of systems 
with an odd number of constituants, providing the un- 



ambiguous evalution of the overlap for the special case of 
two quasiparticle vacua linked by a symmetry transfor- 
mation. 

In the present paper, we present the extension of this 
scheme to the calculation of the overlap between two pos- 
sibly different arbitrary quasiparticle vacua, generalizing 
the treatment of completely filled single-particle states to 
the most general case. This extension thus allows to han- 
dle quasiparticle vacua obtained from blocked 1-, 2-, . . . n 
arbitrary quasiparticle states. To this aim, we present an 
alternative derivation of the overlap that makes use of an 
extension of the standard Thouless parameterization of 
quasiparticle vacua [11] that is advantageous in the pres- 
ence of completely filled single-particle states and allows 
to avoid many of the matrix manipulations elaborated in 
Ref. [§]. Also, our final expression Eq. (63 1 allows for 



the calculation of the overlap of two quasiparticle vacua 
that are expressed in different single-particle bases that 
do not span the same sub-space of the Hilbert space of 
single-particle states, a situation frequently encountered 
in symmetry-restored GCM codes that use a coordinate 
space representation of the quasiparticle vacua in terms 
of their canonical single-particle bases [T2Tfro] . 

The article is organized as follows: Section II intro- 
duces a generalization of the Thouless parameterization 
of quasiparticle vacua for blocked states that will turn 
out to be useful for the purpose of our paper. Sec- 
tion III reviews key properties of Fermionic coherent 
states and Grassmann calculus that will be needed lat- 
eron and thereby introduces our notation. Section IV 
describes the calculation of the overlap, and Section V 
presents some illustrative examples of overlaps from re- 
alistic calculations. Finally, Section VI summarizes our 
findings. An appendix gives the representation of de- 
terminants and Pfaffians of skew-symmetric matrices in 
terms of integrals over Grassmann variables. 



II. PARAMETERIZATION OF 
QUASIPARTICLE VACUA 



1 The possibility to use Pfaffians for this purpose was already 
conjectured much earlier by Balian and Brezin in Ref. [H], but 
never worked out. 



Let {a\a} and {b\b} be two not necessarily equiva- 
lent single-particle bases of dimension n (n even), with 



2 



which are defined two (not normalized) quasiparticle 
vacua \4> a ) and \<pb) through the parameterizations 2 



n ■ 



(1) 



where c is either a or b, \— ) is the bare vacuum of 
single-particle operators, and where the M^ c ' are skew- 
symmetric matrices. Their elements with indices belong- 
ing to the ensemble of fully occupied states {J°} of car- 
dinality 3 = r c are (can be) put to zero. This consti- 
tutes a natural way to regularize the matrix in the 
presence of fully-occupied states. 

The triangle pointing to the right on top of the product 
sign means that it is a "direct" product, 



n 



At at 
-1/ C (5i 



with jj, > ■ ■ ■ > v > 5 , (2) 



as opposed to a "reverse" product obtained, for example, 
by taking the adjoint of Eq. i.e. 

< 

J~[ r, — cgc u • • • c M , with fi> ■ • ■ > u > 5 . (3) 

The parameterization Eq. ([I]) of quasiparticle vacua is 
not very different from the one by Thouless How- 
ever, it has two advantages important for our purpose. 
First, this parameterization is well-defined when dealing 
with fully occupied states. And second, it allows to pa- 
rameterize a quasiparticle vacuum for systems with odd 
particle number in terms of single-particle states, some- 
thing that cannot be achieved with the standard Thou- 
less formula. However, this parameterization can be set- 
up only in a specific single-particle basis that separates 
the fully-occupied single-particle states from the others. 
Such a single-particle basis is, for example, the canonical 
single-particle basis of a quasiparticle vacuum. The use 
of these bases does not impose a serious restriction, since 
they provide the most compact representation of a quasi- 
particle vacuum, such that their use is often desirable in 
numerical applications. 

Finally, the convention for the ordering of single- 
particle levels in Eq. ([!]) (matrix elements of and 
product ordering) is to be kept unchanged for each calcu- 
lation involving a given state \4> c )- In fact, Eq. is just 
an alternative convention that circumvents the conven- 
tion of Ref. [7] to connect the phase of ((j>a\<i>b) to (<j) a \— ) 
and (— \(j>b), which cannot be achieved when either {<j) a \— ) 
or (—\<j>b) (or both) is (are) zero. 



2 For time being, n will be the minimal number of single-particle 
states that allows to represent both \<t> a ) and \<j>b) in their respec- 
tive single-particle basis. 

3 The cardinality is the number of elements of an ensemble. 



III. A REMINDER ON FERMIONIC 
COHERENT STATES AND GRASSMANN 
CALCULUS 

In order to evaluate the overlap (4> a \<l>b) between the 
two states, we introduce, following closely Ref. [7J, two 
sets of Fermionic coherent states 



\z c ) = e" ■"" 
(zc\ = He" 



(4) 
(5) 



for c — a,b, parameterized in terms of anticommuting 
z Ck and z Ck elements of a Grassmann algebra (/, where 
the notations c 1 .z c and z c .c used in Eq. ^ and Eq. ^ 
stand for 



.z r 



i=l 

n 



(6) 
(7) 



In particular, we notice that the coherent states \z c ) are 
not normalized. Instead, one has (— \z c ) = (z c \— ) = 1. 

In what follows, we recall some useful properties 
of Grassmann algebra, its associated calculus, and of 
Fermionic coherent states that will be needed for the for- 
mal derivations outlined below. Concerning Grassmann 
algebra and calculus [TTttTO] . we recall that 

• The adjoint operator performs a one-to-one map- 
ping within Q 



( z c k z ci ) — z c t z c 



Grassmann variables anticommutc 



= 0. 



(8) 
(9) 

(10) 



(11) 
(12) 

(13) 
(14) 



In the following, a product of p Grassmann vari- 
ables (^/-variables) will be called a monomial of de- 
gree p. When p is even (odd) , such a product will be 
called an even (odd) monomial. We notice that an 
even monomial of Q -variables commutes with even 
and odd monomials of ^-variables. We also remark 
that exponentials of pairs of (/-variables also com- 
mute with even and odd monomials of (/-variables, 
such an exponential being a sum of even monomials 
of (/-variables. 

• (/-variables commute with complex numbers and 
anticommute with Fermionic operators. 
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• The fundamental Grassmann calculus rules are 

d 



dz Ck z Ck 



1. 



dz Ck = 0. 



(15) 
(16) 



As a consequence of Eq. (14 1, the Cf-variables z, 



play the role of their own 5-functions, e.g., for an 
analytic function / of (/-variables (c.f. QI5], p. 91), 
we have 



dz Ck z Ck f{z Ck ) = /(0). 



(17) 



• The adjoint variables z Ck and z Ck are independent 
integration variables (c.f. [T7], P- 28). 

Concerning Fermionic coherent states, the properties 
to be used in what follows are: 

• They are eigenstates of second quantized operators 



Cfc \z c 



(Z r . Z c 



(18) 



with the eigenvalues being Grassmann variables. 

• They resolve the identity through the closure rela- 
tion 



t n = I d£z c d™2 c \z c 



(19) 



where l c means that the resolution of the iden- 
tity is built for Fock spaces generated by {c*,c}. 
We use the short-hand notation z c .z c = z\z c — 
J27=i ZciZci- Finally, d" and d™ represent prod- 
ucts of differential elements that are ordered such 
that 



d"a; = dx n ■ ■ ■ dx-i dx\ 
d"a; = dx\ dx2 ■ ■ ■ dx n 



(20) 
(21) 



As compared to [7] and many textbooks, we change the 
ordering of the products of differential elements in the 
integral in order to anticommute them in a more trans- 
parent way, i.e. we use 



d n z d n z 



(22) 



which is equivalent to the more commonly used ordering 4 



(23) 



IV. EVALUATION OF THE OVERLAPS 

A. Preliminary considerations 

To evaluate the expression for the overlap, we start by 
inserting two closure relations, the left (right) one being 
based on the single-particle basis of the left (right) state, 

e.g. {(t> a \4>b) = ((/)a\t a lb\(f)b) 



>,,; = I dgz a dZz a d2z b d%z b 
(za\z b )e 



Zb\<Pb) 



(24) 



where we implicitely use that d™z c d"z c (c = a, b) are 
even (2n) monomials of Grassmann differential elements, 
and thus commuting with Fermionic operators, in order 
to move all differential elements to the very left. 

When the single-particle bases of a and b do not span 
the same subspace of the Hilbert space of single-particle 
states, e.g. when they are not linked through a unitary 
transformation, the resolution of the identity l a (TL b ) 
works for the left (right) state alone, and the two clo- 
sure relations are not equivalent. However, there is no 
loss of generality for the following, their non-equivalence 
being carried by the overlap kernel (z a \z b ). 

Given that by definition (— \z c ) = (z c \— ) = 1, we first 
evaluate the three overlaps 

(4> a \z a ) = ft Zai e-^< >Z ^, (25) 

> 

(z b \ct>b) = e+iZ*' M M J] z b] , (26) 

mi"} 

(z a \z„) = e^M^^H. (27) 

The two first expressions use the properties Eq. ( fi"8| of 
coherent states. The last one uses the Baker-Campbell- 
Hausdorff formula, 5 such that 



{Za\Z b ) 



|e*- 4 e 6t •*»!-) 

, e b i . Zbe 2 a .a e [z a -&,tf. Zb ] 



(28) 



Indeed, the latter is applicable because the commutator 



z a .a,b .z b 



= ^z Q . {ai,St}z b . (29) 

ij 

= y^z ai RijZ bj = z\Rz b , (30) 



commutes with z a .a and I) .z b , where Rij = {a^S^} de- 
notes the matrix of overlaps of the single-particle states 



4 In the commonly used ordering, there is no need to define a 
particular product ordering as dz Ci dz Ci are even monomials of 
5-variables, and the overall order convention is carried only by 
dz Ci dz Ci . 



5 The Baker-Campbell-HausdorfT formula states that, if 
[X, [X,Y]] = [Y, [Y,X]] = 0, then e x e Y = e Y e x e( X ' Y l 
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corresponding to fij and St. We finally obtain Eq. (27) 
by considering that |— ) is a common vacuum for the op- 
erators a and b. 



B. Integration of the reproducing kernel 

As the next step, we integrate the expression for the 
reproducing kernel 

e-*-*"<z„|z 6 )e-*»-*» = e-*-*^^^*-**-**, (31) 

where we use that exponentials of pairs of ^-variables 
commute, thereby allowing to merge the three exponen- 
tial factors. 



Noticing that, in Eq. (24), z a and z b only appear in the 
reproducing kernel Eq. (31), we want to integrate these 



variables separately. In order to do so, we first remark 



that the expression Eq. (31) can be moved to the very 



right of Eq. (24) because it commutes with tj-variables. 



We can as well move the product of differential elements 
d™z a d™z b in front of it by virtue of 



{-l) 2n 4»z a a»z b d»z a d$z b 
d»z a d$z b dSz a d»z b . (32) 



With d"z a d™z b being an even product of Grassmann 
differential elements, it commutes with (<p a \z a ) and 
(z b \(j) b ), quantities containing neither the variables z a nor 
z b , c.f. Eqns. (25) and (26). Rewriting Eq. (24) in a more 



suitable way now gives 



(<t>a\<f>b) = / d3* a d£* 6 



{(j) a \z a )(z b \(l) b ) 



V 



v 



jn - jn -2: a .js a +2;* Rz b -z b .z b 
u > ^>a u <] / 'O e 

v 

reproducing kernel integral 

6,7 
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(33) 



We now evaluate the reproducing kernel integral. Provided that R is non-singular, ' we make the change of variables 
(c.f. Ref. PU p. 14) 



rj f = z\-z\R- x 
t) = z b - R 1 z a 

such that 

ffRrj = -z a .z a + z\Rz b - z b .z b + z t b R~ 1 z a . 
The reproducing kernel can now be written 

g— Z a .Z a +%a RZb-Zb-Zb — p - R~ 1 z a grf Rr] 

The Jacobian of the transformation being one, i.e. d™z a d™z b = d™fjd™r], the integration gives 



J»5 ~z a .z a +z t a B,z h -z b .z b 
^ **a u <i <*o e 



''J d^d^e^^ 

= (-l)"det (R) e - gt > R ~ lz " 
where we have used the determinant formula outlined in Eq. ( A2 ) of Appendix [A] 



(34) 

(35) 
(36) 



(37) 



6 The case of singular R is not equivalent of having zero overlap. 
As an example, consider the case of a partially or completely 
empty single-particle state of the left vacuum which is orthogonal 
to all single-particle states of the right vacuum. In this case, R 
is singular, whereas the overlap is not necessarily zero. 

7 In case of singular R, one has to complete the single-particle 
bases of \<f> a ) and \cj> b ) in order to get an invertible matrix, for 
example using a Gram-Schmidt orthonormalization procedure. 
Still, having a non-singular matrix R is not equivalent to having 
equivalent bases a and b. 
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C. Re-expression of the overlap 



Using Eqns. (25), (26) and (37), we are now able to rewrite Eq. (33) as 



(4> a \<k) = (-l)"det(i?) / dX^t^NaX^l^e" 1 ^) (38) 
= (-l)"det(i?) / dXCzi ft *«< fl z- bj e(^ E -< , * z - 2 - + ^-< ) «-^)) (39) 

where we used that exponentials of pairs of (/-variables commute with (/-variables. Closely following the notation of 
Ref. [7J, we introduce the matrix M and the vector £ 



such that 



1 M 

2 



(40) 

(41) 
(42) 



Inserting these definitions into Eq. ( 39 ) , the overlap kernel becomes 



(-l) n det(R) f d3* d£« 6 ( f[ z CH fl 



Zb 4 



C*mc 



(43) 



D. Integration over fully-occupied states 



We will now integrate variables corresponding to fully occupied states in Eq. (43) by virtue of Eq. (17). In order 



to do so, we need to move the corresponding differential elements to the very right, which gives a sign factor because 
of the anticommutation of Grassmann differential elements. For example, such a rearrangement for a single variable 
in d™z c and d™z c gives 



d"z c = dz Cl ■ ■ ■ dz Ck ■■ - dz, 



1 



Udz Cz ] (-ir- k dz Ck , 



(44) 



d > z c = dz Cn ■ ■ ■ dz Ck ■ ■ ■ dz Cl 



]Jdz Ci ) {-if-Uz. 



When there is more than one fully occupied state, the repeated application of this procedure gives 



(45) 



dx = ( n dz c% n {~i) n - k dz Ck , 



d^Zc 



n d ^ n i- i ) k ~ i ^c 

iW} J \ke{i°} 



(46) 
(47) 



6 



where in each equation we notice the opposite order for the products over indices corresponding to fully occupied 



states as compared to other states. Combining Eqns. (46 1 and (471 gives 

X — 



d> a d^ b = [ n dz «A n {-^ n - k dza k \ n ^ n 

mi a } J \k£{i-} J W{/ 6 } J \ie{i»} 



(48) 



= (-i) nr ° ( n **« n d ^ n (-v 1 - 1 ^ n (-ir - **** i • <- w > 

mi*} J \mi"} J \mi"} J Wun 



The additional sign in Eq. (49 1 comes from the commutation of variables k £ {I a } to the very right, as indicated in 

(50) 



Eq. ( 48 ) . Defining the sign factor 



Eq. (49) finally gives, after a suitable rearrangement of the sign factors 



d> d^6=(-ir r " +r ^[ n n dz »A n ^ n 

mi*} J XJW} J Veil"} J \ke{i«} 



(51) 



The ordering of the differential elements corresponding to indices of fully occupied states are now in the appropriate 
order with respect to their associated products in Eq. (43) to perform their integration. 



By virtue of Eq. (17 1, the integration over fully occupied levels of the integrand in Eq. (43) gives 



/ ( n *.) ( ri *») ( ri *. n *» #«) - (^-'U^,,,,,^,, 

\ke{i b } ) \fce{/-} / \ie{i«} je{i"} J x *' x ' 



(52) 



With this, Eq. (43) can be rewritten as 



(^|^) = (-l)"(-ir» +r Vdet(i?) / I n dz a ) I f[ dz b \ (e^) (53) 

J \mi*} J \mh J K h^mi>hn-mi*} 

where the integration only runs over variables of indices associated to not fully-occupied states. 

E. Integration over the remaining variables 

We now define a new matrix M r , sub-matrix of M where rows and columns of indices i £ {I b } of fully occupied states 
in \(f>b) and of incices j + n, where j € {I a } are fully occupied levels in \<p a ), have been removed. The corresponding 
appropriate vector £ r is built from £ in the same manner, removing components with indices i £ {I b } and j + n such 
that j £ {I a }- The matrix M r is skew-symmetric with dimension iV r x N r , and the vector £ r has N r elements, with 
N r = 2n- (r a + r b ). 

For an example where there are two fully occupied states in each quasiparticle vacuum, with indices i, k for the 
state \(f>b) and indices j, I for the state \4>a), respectively, the matrix M r and the vector £ r can be schematically 
represented as 



ik 



J 1 









\ 


-f- 






-R- 1 










V 


(« ')' j 




(MW)* 

) 



Cr - j 



Zb 



Z a 

V J 



(54) 
(55) 



(56) 
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where labeled rows and columns have been removed from the original objects M and £. Using the submatrix M r and 
subvector £ r , Eq. (531 can be rewritten as 



b a \<t> b ) = (-l)"(-l)" r " +r Vdet(i?) 



/ f ri d Za ) ( n *J 



(57) 



In order to apply the Pfaffian formula Eq. (A4), the order of n^f/"} dz ai has to be reversed, such that the differential 
elements in Eq. (57 1 are in the appropriate order with respect to the matrix M. r . This is achieved by reversing the 
order of dz ^ 



n dz ai = j| dZai . 



(58) 



The sign factor from the reversal of the product Z\ ■ ■ ■ z n = (— x )/ 2 z n ■ ■ ■ Z\ can be obtained by induction. The 
differential elements are now in the appropriate order 



such that we can now integrate the remaining variables using the Pfaffian formula, Eq. ( A4 1 , 

((t> a \<f> b ) = (-l) n (-l) nr « +r »a dot (R) (_i)(»-'--)("-'---i)/2 J d^C r e^' M ^ 

= (_l)»('H-l)/2(_ 1 )r.(r.-l)/aj_ 1 jr a +r fcCTdet ^ pf ^ Mr j < 



(59) 
(60) 



(61) 
(62) 



This formula, however, assumes N r to be even, see ap- 
pendix |Aj When N r is odd, \4> a ) and \4>b) have in fact 
different number parity [20] , In that case, the overlap 
(0a 1 0b) is automatically zero. From the definition of the 
Pfaffian, which by definition is zero for skew-symmetric 
matr ices of odd rank, we can thus notice that formula 
Eq. (57) can still be applied. In particular, (— l) r «+ ri > 
will always be one except when multiplied by zero, al- 
lowing to drop this sign factor in the final formula. We 
thus summarize the final expression for the overlap as 



where 



(<M<M = s„s ra o- det (R) pf (M r ) 



Sra = (-ly.c--- 1 )/ 2 , 

a — <_ ^Efe=l..r a ( n +«o)+Efc=l. 



(63) 



(64) 
(65) 
(66) 



The sign factors depend on the number of states in 
the single-particle bases n, the number of fully occupied 
states r a in a, and the indices ik of fully occupied states 
in the bases c — a, b. 



Equation (63) provides the generalization of Eq. (7) 



vacua with an arbitrary number of fully occupied single- 
particle states. In particular, it can be applied to over- 
laps that involve an odd number of blocked quasiparticle 
states, a case not considered at all in Refs. [71 [9], and for 
the special case of symmetry restoration only in Ref. [TU] . 
Moreover, when an even number of particles is fully oc- 
cupied (either for blocked 2n quasiparticle states, or as a 
result of the minimization, or both), Eq. (63) provides a 



formally justified alternative to the regularization of the 
matrix M performed in [9] . 

Besides this important generalization, there is another 
noteworthy difference to previous work by Robledo [7J [S] . 
Indeed, Eq. ( 63 ) is directly expressed in the single- 
particle bases of \<f> a ) and \4>b), respectively, that allow 
for the most compact representation of these quasipar- 
ticle vacua. In particular, Eq. ( 63 ) can also be applied 



of Ref. [7 to the overlap between different quasiparticle 



without invoking a complete single-particle basis span- 
ning the single particle subspace a U b, cf. the discussions 
above. However, as explained there, if the matrix R is 
singular, one is forced to complete each single-particle 
basis until det (R) is achieved. 

For quasiparticle vacua for which there are no fully 
occupied states in their respective single-particle basis, 
it is easy to show that Eq. (7) of Ref. [7j is recovered. 
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o.o 

Re ((0 a 1 0b)) 

FIG. 1: Real and imaginary parts of the overlap without 
particle-number projection for the lowest one-quasiparticle 
state in 25 Mg obtained with SHI. The Euler angles a and 
P are held fixed at values of a — 1.25° and /3 = 7.17°, wheras 
7 is varied in the interval [0, 720°] with a discretization of 
288 points. Filled circles on the curve represent a discretiza- 
tion of 48 points in the interval [0, 720°], which is sufficient to 
converge observables. Note the difference in scale of real and 
imaginary parts. 
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FIG. 2: Same as Fig.[T] but for a = 43.75° and /3 = 71.94°. 



V. SOME ILLUSTRATIVE EXAMPLES 

The determination of the phase of the overlap by the 
widely used techniques that rely on a Taylor expansion 
of the overlap around a matrix element of known phase 
[6l 114] works very well when restricting the calcula- 
tions to time-reversal invariant HFB states. However, it 
becomes increasingly difficult when time-reversal is bro- 
ken, as it happens for odd-^4 or odd-odd nuclei, or for 
states obtained with cranked HFB, cf. the examples dis- 
cussed in Refs. P2H [21]. 

We have implemented Eq. ( 63 ) into our numerical 



codes for particle-number and angular-momentum re- 
stored GCM calculations based on triaxial HFB states 
[THUS], using the routines for the calculation of the Pfaf- 



fian of Ref. [22J . We now present three examples where 
techniques to follow the overlap through Taylor expan- 
sion might fail and the direct calculation of the overlap 
becomes a necessity. These illustrations will show trajec- 
tories in the complex plane of overlaps of quasi-particle 
vacua as obtained during angular-momentum projection 



(0 a |0 b ) = (<f\R(a,l3,j) \cp), 



(67) 



where R is the rotation operator that depends on the 
three Euler angles a, f3, and 7. 

The first two examples are presented in Figs, [l] and [2] 
They illustrate the trajectory of the overlap in the com- 
plex plane when varying the Euler angles 7 when pro- 
jecting the lowest self-consistent one-quasiparticle state 
of 25 Mg, for two different combinations of a and j3. The 
first one, Fig. [I] illustrates that real and imaginary parts 
of the overlap can vary on quite different scales. In this 
particular case, most of the modulus of the overlap is 
carried by the real part, and the phase of the overlap 
is most of the time either close to zero or close to ±7r. 
Unless the discretization of the Euler angles is carefully 
adapted, the phase of the overlap might change by al- 
most 7r when crossing the imaginary axis, which is very 
difficult to distinguish from a discontinuity by 7r encoun- 
tered when having lost the phase. The second example, 
Fig. [2j obtained for a different combination of Euler an- 
gles a and (3, shows that the trajectory of the overlap in 
the complex plane may exhibit cusps, which might again 
be difficult to resolve when discretizing Euler angles. 

In Fig. [3J we present the trajectory of the overlap 
of a high-spin state in 24 Mg, obtained from cranked 
HFB+Lipkin Nogami [23] • The two inserts illustrate that 
variations may occur on very different scales with quite 
involved structures. 

These three examples demonstrate that a Taylor- 
expansion-based algorithm to determine unambiguously 
the sign of the overlap may become difficult in applica- 
tions that use time-reversal invariance breaking quasi- 
particle vacua. Indeed, it should be able to resolve dis- 
continuities or cusps, or other involved structures. The 
latter might happen at very different scales, and the dis- 
cretization must be chosen in order to account for all 
these details. Furthermore, we expect the complexity 
of such trajectories to increase with increasing intrinsic 
angular momentum. A direct determination of the over- 
lap is thus a considerable improvement not only from a 
formal point of view, but also from the perspective of 
the complexity of reliable algorithms for the computa- 
tion of the overlap on the one hand and of computing 
time on the other hand, as it will often allow for the use 
of a smaller number of combinations of Euler angles in 
angular-momentum projection. 



VI. DISCUSSION AND OUTLOOK 



To summarize our main findings 
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FIG. 3: Real and imaginary parts of the overlap without 
particle-number projection for the cranked I = 8h state in 
24 Mg obtained with SHI. The Euler angles a and /3 are held 
fixed at values of a = 23.75° and /3 — 66.96°, respectively, 
wheras 7 is varied in the interval [0, 360°] with a discretiza- 
tion of 144 points. The inserts amplify the zone at very small 
overlaps, whereas filled circles on the curves represent a dis- 
cretization of 24 points in the interval [0, 360°], which is suf- 
ficient to converge observables. 



1. An extension of the Thouless parameterization of 
quasiparticle vacua with completely filled single- 
particle states allows to calculate the overlap di- 
rectly in a formalism based on Grassmann algebra 
and coherent states similar to the one outlined in 
Ref. [7] , but in such a manner that fully occupied or 
empty single-particle levels are automatically taken 
care of without any need for the manipulation (or 
regularization) of matrices elaborated in Ref. [5] . 

2. This extension of the Thouless expression allows 
to handle all possible quasiparticle vacua that 
have completely filled states, i.e. also 1-, 2-, ... n- 
quasiparticle states, not just quasiparticle vacua 
that can be expressed as limits of fully-paired quasi- 
particle vacua as in Ref. [9]. 

3. The handling of blocked states is not restricted to 
pure symmetry restoration as the one proposed in 
Ref. [10], and therefore can be also applied when 
the non-rotated left and right states are different, 
which is necessary for GCM calculations. 

4. Our final expression for the overlap allows for the 
calculation of the overlap of two quasiparticle vacua 
that are expressed in two different single-particle 
bases that do not span the same sub-space of the 
Hilbert space of single-particle states. The knowl- 
edge of a complete basis spanning both single- 
particle bases is not needed, as compared to Ref. [S]. 
In this way, the technique can be directly imple- 
mented in codes that use a coordinate space rep- 
resentation of the quasiparticle vacua in terms of 
their canonical single-particle bases [T3TU6) . 



The expression has been implemented into our numer- 
ical codes for particle-number and angular-momentum 
restored GCM calculations based on triaxial HFB states 
using the full space of occupied single-particle states 
[TBI ITT]] . It has been extensively tested for symmetry 
restoration and for the calculation of non-diagonal ma- 
trix elements for symmetry-restored GCM calculations 
without encountering cases where it fails. 

By contrast, the technique "to follow the phase" from 
Ref. [S] H3] , or the one to follow the overlap in the com- 
plex plane from Ref. [5] require often to use a very fine 
discretization to resolve the sign ambiguity when pro- 
jecting on angular momentum as soon as time-reversal 
invariance of the HFB states is broken. In particular, 
it may become necessary to use a discretization of the 
integrals over Euler or gauge angles that are much finer 
than what is actually needed to converge observables. In 
addition, the direct calculation of the overlap also has 
the advantage to avoid complicated coding for the set-up 
of a reliable Taylor-expansion-based algorithm, in par- 
ticular since reliable routines to compute pfaffians are 
available [55]. Finally, in the general case of configura- 
tion mixing where there might not be a symmetry that 
establishes a reference sign for the overlap, such a direct 
calculation of the overlap could become to be mandatory. 

In summary, we report an expression for the overlap 
between arbitrary quasiparticle vacua that is easy to cal- 
culate and that is very robust in realistic applications. It 
is a key ingredient for the extension of symmetry restora- 
tion and Generator Coordinate Method-type calculations 
to angular-momentum-optimized states, either by crank- 
ing, or by blocking. 



Acknowledgments 



This work was supported by the Agence Nationale de 
la Recherche under Grant No. ANR 2010 BLANC 0407 
" NESQ" . The computations were performed using HPC 
resources from GENCI-IDRIS (Grant 2011-050707). 



Appendix A: Some remarkable Gaussian integrals 
over Grassmann variables 



For the readers' convenience, we give the identities that 
represent determinants and pfaffians as integrals over 
Grassmann variables using our convention for the order 
of the differential elements, Eq. (22). 
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1. Determinant 

The first one is the determinant identity that is defined, 
for a given matrix M, as 

det(M) = [ J| (dz.dzi) exp ( ^ ZiM ijZj \ (Al) 

= (-1)" J d£zd>exp (f^ Zi M i:i z)jh.:i) 

where the first equation is the expression from Ref. [19] , 
p. 13, Eq. (1.67), and the second uses an alternative con- 
vention in the ordering of differential elements. The lat- 
ter differs by a sign because of the anticommutation of 
Q- variables. 



2. Pfafflan 

The second one is the Pfafhan identity that is defined 
for a skew-symmetric matrix A of dimension 2n x 2n, as 



(see [IS], p.15, Eq. (1.80)): 



pf (A) = J dz 2n ■ ■ ■ dz 2 dz x exp f i ^ ZiAijz\k$) 
(l^iAijzA ■ (A4) 



d^™2;exp 



To obtain the correct sign, it is crucial that the differen- 
tial elements are in the same order as the indices of the 
matrix Ay. This implies sometimes to make some ma- 
nipulations to bring the entire expression into the proper 
form, as for example between Eq. (57) and Eq. (61). 



As can be easily seen, this definition is only valid for 
matrices A of even rank. The determinant and the Pfaf- 
fian of a skew-symmetric matrix A are related by 



[pf(^)] 2 = det (A) 



(A5) 



As the determinant of a skew-symmetric matrix of odd 
rank is always zero, the pfaffian of such a matrix is de- 
fined to be zero as well. 
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